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ABSTRACT 

The dynamical equations describing the evolution of a self-gravitating fluid can be 
rewritten in the form of a Schrodinger equation coupled to a Poisson equation deter- 
mining the gravitational potential. This approach has a number of interesting features, 
many of which were pointed out in a seminal paper by Widrow & Kaiser (1993). In 
particular we show that this approach yields an elegant reformulation of an idea due 
to Jones (1999) concerning the origin of lognormal intermittency in the galaxy distri- 
bution. 

Key words: Cosmology: theory ~ galaxies: clustering ~ large-scale structure of the 
Universe. 



1 INTRODUCTION 

Galaxy redshift surveys such as the 2dF Galaxy redshift sur- 

■ vey and the Sloan Digital Sky Survey are poised to yield a 
copious harvest of statistical information about the distri- 
bution and dynamical properties of the large-scale structure 

■ of the Universe. At the same time relatively few statistical 
properties of the structure revealed by these observational 

, programs are understood in quantitative detail using ana- 

■ lytic methods. 

One particular facet that has received some attention 

■ over the years has been the property that the one-point prob- 
' ability distribution of density fluctuations p{p) appears to 

have a roughly lognormal form, i.e. log p has a roughly nor- 
mal distribution (Coles & Jones 1991). It is now established 
that this property has an interesting connection with the 
scaling properties of moments of the probability of the dis- 
tribution. Taking a generic random variable X, such that the 
distribution of X within cells of side L is denoted p{X; L), 
then the g-th moment at a given value of L is said to display 
scaling if 

=^p(X;L)X'' cxL^''). (1) 

This means that different powers q of the density field vary 
as a different power of the coarse-graining scale L. The func- 
tion iJL{q) is called the intermittency exponent, and it can be 
extracted from observations. Jones, Coles & Martinez (1992) 
showed that observations suggest a roughly quadratic depen- 
dence of pL{q) upon q and that this is related to the under- 
lying near-lognormal form of the density fluctuations. A set 
displaying the form (1) is usually termed a multifractal; see 
Paladin & Vulpiani (1987) for general discussion. 



The term intermittency was coined to described proper- 
ties of stochastic processes described by high skewed proba- 
bility distributions with very slowly decaying tails. The log- 
normal is a prime example, used in a pioneering paper by 
Kolmogorov (1962). In the time domain it refers to pro- 
cesses involving extended quiescent periods interrupted by 
bursts of intense activity. In the spatial domain, intermittent 
processes are ones in which isolated regions of high density 
are separated by large voids. Although in a qualitative sense 
the appUcation of the concept of intermittency to large-scale 
structure seems plausible, a quantitative description of how 
it arises is not easy to obtain. Drawing on ideas discussed 
by Zel'dovich et al. (1985, 1987), Jones (1999) suggested an 
analytical model for the cosmological context. It is the pur- 
pose of this paper to describe an alternative formulation of 
the Jones (1999) model and substantially strengthens the 
physical understanding of this model. This also provides an 
opportunity to advocate the wider use of an alternative for- 
mulation of the gravitational instability scenario discussed 
by Widrow & Kaiser (1993). 



2 THE FLUID APPROACH TO STRUCTURE 
FORMATION 

2.1 Basics 

In order to understand how the intermittent form of large- 
scale structure arises, it is best to begin with the stan- 
dard fluid-based approach to structure growth. In the stan- 
dard treatment of the Jeans Instability one begins with 
the dynamical equations governing the behaviour of a self- 
gravitating fluid. These are: the Euler equation 
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(2) 



d(\r) 1 
^ + (v • V)v + - Vp + = ; 
ot p 

the continuity equation 

g + +V(pv) = (3) 

expressing the conservation of matter; and the Poisson equa- 
tion 

\7^(l> = 4ttGp, (4) 
describing Newtonian gravity. 

2.2 The Cosmological Setting 

If the length scale of the perturbations is smaller than the ef- 
fective cosmological horizon (Ih = c/Ho, a Newtonian treat- 
ment of cosmic structure formation is still expected to be 
valid in expanding world models. In an expanding cosmo- 
logical background, the Newtonian equations governing the 
motion of gravitating particles can be written in terms of 



r/ait) 



(5) 



(the comoving spatial coordinate, which is fixed for observers 
moving with the Hubble expansion), 

V = r — Hr — ax (6) 

(the peculiar velocity field, representing departures of the 
matter motion from pure Hubble expansion), p(x,t) (the 
matter density) , and (^(x, t) (the peculiar Newtonian gravi- 
tational potential, i.e. the fluctuations in potential with re- 
spect to the homogeneous background) determined by the 
Poisson equation in the form 



Vx^<?f> = i-JvGa^ip - po) = A-JvGa^poS. 



(7) 



In this equation and the following the suffix on Vx indicates 
derivatives with respect to the new comoving coordinates. 
Here po is the mean background density, and 



5 = 



Po 



po 



(8) 



is the density contrast. Using these variables the Euler equa- 
tion becomes 



9(av) 



+ (V ■ Vx)v : 



P 



(9) 



The first term on the right-hand-side of equation (^) arises 
from pressure gradients, and is neglected in dust-dominated 
cosmologies. Pressure effects may nevertheless be important 
in the the (collisional) baryonic component of the mass dis- 
tribution when nonlinear structures eventually form. The 
second term on the right-hand side of equation (^ is the 
peculiar gravitational force, which can be written in terms 
of g = — Vx0/a, the peculiar gravitational acceleration of 
the fluid element. If the velocity flow is irrotational, v can 
be rewritten in terms of a velocity potential (j>v'- 



-Vx0i>/a 



(10) 



This is expected to be the case in the cosmological setting 
because (a) there are no sources of vorticity in these equa- 
tions and (b) vortical perturbation modes decay with the 
expansion. Next we have the revised continuity equation: 



^ + 3Hp+-V^{pv) 
ot a 



0, 



(11) 



2.3 The Nonlinear regime 

The strength of the fluid approach is that it furnishes a rel- 
atively simple perturbative treatment for the early stages 
of structure formation, wherein the density fluctuations are 
small. Perturbation theory fails when nonlinearities develop, 
and unfortunately this it is in the nonlinear regime where 
scaling and intermittency arise. It is important to stress, 
however, that the fluid treatment is intrinsically approxi- 
mate anyway. A fuller treatment of the problem requires a 
solution of the Boltzmann equation for the full phase-space 
distribution of the system /(x,v,t) coupled to the Poisson 
equation (^ that determines the gravitational potential. In 
cases where the matter component is collisionless, the Boltz- 
mann equation takes the form of a Vlasov equation: 



di 

dt 



3 



^ \ dxi dvi 



' dxiJ 



(12) 



The fluid approach can only describe cold material where 
the velocity dispersion of particles is negligible. But even if 
the dark matter is cold, there may be hot components of 
baryonic material whose behaviour needs also to be under- 
stood. Moreover, the fluid approach assumes the existence 
of a single fluid velocity at every spatial position. It there- 
fore fails when orbits cross and multi-streaming generates 
a range of particle velocities through a given point. While 
the existence and location of the caustics that occur and 
shell crossing can be determined using the Zel'dovich (1970) 
approximation (Coles, Melott & Shandarin 1993), the dy- 
namics and distribution of matter after shell-crossing is not 
described in this approach. Attempts to understand proper- 
ties of non-linear structure using the fluid model therefore 
resort to further approximations (Sahni & Coles 1995) to ex- 
tend the approach beyond shell-crossing. Even the numerical 
N-body methods that many regard as the ultimate approach 
to this kind of problem are also approximate, corresponding 
to a particular type of Monte Carlo integration of the Vlasov 
equation. 



2.4 The Jones Model 

It is within the overall framework of the fluid model that 
Jones (1999) sought to understand the observed intermit- 
tency of the large-scale structure of the Universe. Using the 
velocity potential introduced above, he first introduces an 
effective Bernoulli equation for the flow: 

90„ (V0.)^ 

where is the actual gravitational potential. This equation 
neglects terms involving pressure gradients as mentioned 
above. To cope with shell-crossing events, Jones (1999) in- 
troduces an artificial viscosity v by adding a term to the 
right-hand-side of this equation: 



d(l>v 
~dt 



2a2 



(14) 



The viscosity is introduced to prevent the particle flow 
from entering the multi-stream region by causing parti- 
cles to stick together at shell-crossing. This is also used in 
an approach called the adhesion approximation (Gurbatov, 
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Saichev & Shandarin 1989). Using the Hopf-Cole transfor- 
mation (jii, = ~2v\ogip and defining a scaled gravitational 
potential via = 2ve we can write the Bernoulli equation 
as 



— = !/V (^ + e(x)^. 



(15) 



This is called the random heat equation, because of the 
existence of the spatially-fluctuating potential term e(x). 
The gravitational potential changes very slowly even during 
nonlinear evolution (Brainerd, Scherrer & Villumsen 1993; 
Bagla & Padmanabhan 1994), so Jones (1999) assumes that 
it can be taken as constant and to be Gaussian distributed. 
An approximate scaling solution to ( [l5| ) can then be found 
using a path integral adapted from that normally used in 
quantum physics (Feynman & Hibbs 1965); see below for 
more details. In this approximation, the function tp has a 
lognormal distribution (Coles & Jones 1991). We refer the 
reader to Jones (1999) for details; see also Zel'dovich et al. 
(1985, 1987). 

This model is one of the few attempts that have been 
made to understand the non-linear behaviour of the mat- 
ter distribution using analytic methods. Although not rigor- 
ous it surely captures the essential factors involved. It does, 
however, suffer from a number of shortcomings. First, the 
approach does not follow material beyond the shell-crossing 
stage. Second, the viscosity v that is needed does not have 
properties that are very realistic physically: it can depend 
neither on the density p nor the position x. Moreover, in 
the final analysis Jones (1999) takes the limit — > 0, so it 
cancels out anyway. One is tempted to speculate that its in- 
troduction may be unnecessary Third, the function (p(x, t) 
that emerges from equation (hsl) is not the desired density 
p(x, t) nor does it bear a straightforward relation to the den- 
sity. Finally, it is not clear how the motion of a coUisional 
baryonic component can be modelled within this framework. 



3 WAVE MECHANICS AND STRUCTURE 
FORMATION 

A novel approach to the study of collisionless matter was 
suggested by Widrow & Kaiser (1993). They advocate the 
replacement of the Euler and continuity equations by a ver- 
sion of the Schrodinger equation 



ih— = - — V ip + m0(x)-(/'. 
at 2m 



(16) 



The equivalence between this and the fluid approaches has 
been known for some time; see Spiegel (1980) for historical 
comments. Originally the interest was to find a fluid inter- 
pretation of quantum mechanical effects, but in this context 
we shall use it to describe an entirely classical system. In this 
spirit, the constant h is taken to be an adjustable parameter 
that controls the spatial resolution A through a de Broglie 
relation A = fi/mv. To accommodate gravity we need to 
couple equation (|l6|) to the Poisson equation in the form 



(17) 



It clearly follows that p — IV'P; the wavefunction ip{^,i) 
evidently complex. Moreover, if the matter distribution is 
initially cold then one can further construct a wavefunction 



that also encodes the velocity part of phase space in its 
argument through the ansatz 



V'(x) = \/p(x)exp[i6l(x/ft)] 



(18) 



where V6'(x) = p(x), the local 'momentum field'. This for- 
malism thus yields an elegant description of both the density 
and velocity fields in a single function. 

This formulation provides a useful complementary ap- 
proach to many techniques, including A^-body methods. It 
also provides a new light with which to study the Jones 
(1999) model. Widrow & Kaiser (1993) show using theoret- 
ical arguments and numerical simulations that this system 
allows accurate numerical evolution of the system beyond 
shell-crossing, so it does not have the ad hoc construction 
needed by the Jones (1999) model to remedy this. 

Second, no artificial viscosity is required. Equations ( p^ 
and (^ are of the same form, apart from minor subtleties 
like the use of complex time coordinates. The potential term 
is easily understood in (|l^) , and the wavefunction tp now has 
a straightforward relationship to p. The upshot of this is that 
if one adopts the approximation of constant gravitational 
potential one can use the same path integral approach as 
described by Jones (1999). In a nutshell, given some initial 
wavefunction ip{x',t') one can determine the wavefunction 
at a subsequent time i^in, t) using 



7/;(x,t)= / K{^,t;x.',t')'4>{x.',t')d'^x', 



(19) 



where the function A'(x, t; x', t') involves a sum over all 
paths r connecting the initial and final states with t > t': 



K{x,t-x.',t') = vrexp[is{r)/h], 



(20) 



where D is an appropriate measure on the set of classical 
space-time trajectories. For a particle moving in a potential 
V{x,t) — m(j){-x.,t) the action S for a given path F is given 
by 



1 .2 



(x) 



dt. 



(21) 



Note the presence of the Gaussian field in equation (21) and 
hence in the exponential of the integrand on the right-hand- 
side of equation (|o|) . To get an approximate solution to this 
system we can follow the same reasoning as Zeldovich et al. 
(1985, 1987) and Jones (1999), ignoring time-varying terms, 
using the Gaussian properties and counting the dominant 
contributions to the path integral to deduce that the inte- 
gral produces a solution of lognormal form. This part of the 
argument is identical to that advanced by Jones, except that 
the solution is for i/) rather than ip and since p is then 
one directly obtains a lognormal distribution for the desired 
density p(x, t). 

It should be stressed that, although the present ap- 
proach clearly provides a more elegant formulation of the 
problem, the deduction of lognormality remains approxi- 
mate; the lognormal is not the exact solution to the system 
to either Jones' equation (15) or the present equation (16). 
How accurately this approximate form applies is open to 
doubt and will have to be checked by full numerical solu- 
tions. Interestingly, however, it is known to apply quite ac- 
curately in quantum systems such as disordered mesoscopic 
electron configurations (Janssen 1998). As mentioned above, 
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the Schrodinger approach yields a wavefunction tp which is 
directly related to the particle density p via p — \ip\^ ■ Such 
quantum systems also display lognormal scaling for proper- 
ties such as the conductance, which depends on It is 
a property of the lognormal distribution that if a random 
variable X is lognormal, then so is X" . In such systems the 
role of the gravitational potential is played by a potential 
that describes the disorder of a solid, perhaps caused by the 
presence of defects. Such systems display locahsation at low 
temperature which is similar in some ways to the original 
idea of Anderson location (Anderson 1958). The formation 
of strongly non-linear structures by gravity is thus directly 
analogous to the generation of localised wavefunctions in 
condensed matter systems. 

Finally, and perhaps most promisingly for future work, 
the equation ( |l^ ) offers a relatively straightforward way of 
modelling the behaviour of coUisional material. The addi- 
tion to the potential of a term of the form ajV'P (with a an 
appropriately-chosen constant), converts the original equa- 
tion (hd) into a nonlinear Schrodinger equation: 

ih^ = —^-V^^ + m4>{:>i)^ + a\-tpf^P (22) 

(Sulem & Sulem 1999). This equation is now equivalent to 
those that describe the flow of a barotropic fluid; see Spiegel 
(1980). This system can therefore be used to model pressure 
effects, which are otherwise only handled effectively using 
numerical methods such as smoothed-particle hydrodynam- 
ical approximations (e.g. Monaghan 1992). In the context 
of quantum systems, the nonlinear term is used to describe 
the formation of Bose-Einstein condensates (e.g. Choi & Niu 
1999 and reference therein). 



4 DISCUSSION 

In this short paper I have sketched out an approach to the 
study evolving cosmological density fluctuations that relies 
on a transformation of the Vlasov-Poisson system into a 
Schrodinger-Poisson system. The transformation is not a 
new idea, but despite the efforts of Widrow & Kaiser (1993) 
it does not seem to be well known in the astronomical com- 
munity. The immediate advantage of this new formalism is 
that it yields a rather more convincing approach to under- 
standing the origin of spatial intermittency and approximate 
lognormality in the galaxy distribution than that offered by 
Jones (1999). It also makes a connection in the underlying 
physics with other systems that display similar phenomena. 

On the other hand, one must be aware of the approxima- 
tions also inherent in the present approach. The Schrodinger 
equation is not exact, and its usefulness as an approximate 
tool is restricted by a number of conditions outlined by 
Widrow & Kaiser (1993); see also Spiegel (1980). Further- 
more, the lognormal solution of the system is a further ap- 
proximation and may not be valid especially in the strongly- 
fluctuating limit. Although it neatly bypasses some of the 
problems inherent in the Jones (1999) model, the nonlinear 
wave equation is by no means easy to solve in general sit- 
uations. Numerical methods will still have to be employed 
to understanding other aspects of the evolution of cosmic 
structure within this framework as indeed they are in other 
branches of physics. 



One particular issue worth exploring using this ap- 
proach is to understand the limits of the approach in 
strongly non-linear situations. As it stands, the justification 
for the lognormal approximation arises from the weakly non- 
linear behaviour of coUisionless matter moving in an almost 
constant potential field. Taking into account the expansion 
of the Universe, the changing gravitational potential, and 
the possible effects of matter pressure within in the action 
formalism described in Section 3 may well reveal that a dif- 
ferent form of hierarchical scaling pertains in the strongly 
non-linear regime. One aspect of this is that the hierarchy of 
correlation functions that describe a lognormal distribution 
display Kirkwood (1935) scaling, while it appears from nu- 
merical A''-body simulations that cosmological fluctuations 
display a different hierarchical form. For a discussion of the 
relationship between lognormal and hierarchical scaling, see 
Coles & Frenk (1991). 

Apart from the particular topic of intermittency, this 
paper is also intended to stimulate interest in this general 
approach. Many techniques exist for studying the wave me- 
chanics of disordered systems, such as the renormalization 
group and path-integral methods, few of which are used by 
cosmologists working in this area. It is to be hoped that 
the introduction of some of these methods may allow better 
physical insights into the behaviour of non-linear structure 
formation than can be found using brute-force A'^-body tech- 
niques. 
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